DESY 04-181 ISSN 0418-9833 

MPP-2004-116 

'hep-ph/0409303] 

September 2004 

Next-to-leading-order predictions for D*^ plus jet 
photoproduction at DESY HERA 

"t;!]-" ■ Gudrun Heinrich 

O . 

O ■ • Institut fiir Theoretische Physik, Universitdt Hamburg, 

(N 

Oh' 
(D 

l> 
(N 



X 



Lumper Chaussee 149, 22761 Hamburg, Germany 

Bernd A. Knieh|*l 

Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut) , 



^ ■ Fdhringer Ring 6, 80805 Munich, Germany 

cn ■ 

^ : Abstract 

' We study the photoproduction of a D*^ meson in association with a hadron jet at next-to-leading 

■ order in the parton model of QCD with non-perturbative fragmentation functions extracted from 



LEPl data of e^e annihilation. The transverse-momentum and rapidity distributions recently 



Q^' measured at DESY HERA in various kinematic ranges nicely agree with our theoretical predictions. 

' This provides a useful test of the universality and the scaling violations of the fragmentation 



functions predicted by the factorization theorem. These comparisons also illustrate the significance 



^ \ of the charm component in the resolved photon. This is elaborated by investigating the cross-section 

distributions in x^-^^ and cos 6**. 
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I. INTRODUCTION 



Heavy-flavor production has always been an important testing ground for quantum chro- 
modynamics (QCD), one of the reasons being that it addresses the problem of where to draw 
the dividing line between perturbative and non-perturbative aspects. The photoproduction 
of open charm is particularly interesting because it also allows valuable insights into the 

axtonic structure of the photon. Results on inclusive D*^ photoproduction from the HI 
and ZEUS 12, yl Collaborations at the DESY ep collider HERA have been compared to 
nnn 

leading-order (LO) and next-to-leading order (NLO) calculations |5|, |6|. More recently, 
data on the photoproduction of a D*^ meson in association with a hadron jet have also 
become available 

Concerning the theoretical treatment of open heavy-flavor production, several approaches 
have been followed in the literature. The QCD-improved parton model implemented in the 
modifled minimal-subtraction (MS) renormalization and factorization scheme and endowed 
with non-perturbative fragmentation functions (FFs), which proved itself so convincingly 



for light-hadron inclusive production 111, also provides an ideal theoretical framework for a 

n n 

coherent global analysis of D- p] and _B-meson data provided that /i ^ niq, where /i is 
the energy scale characteristic for the respective production process and Q = c,b. Then, at 
LO (NLO), the dominant logarithmic terms, of the form a" In" {^fj? /mo) (a"'*'^ In" (/i^/mg)) 
with n = 1,2, . . ., where is the strong-coupling constant, are properly resummed to all 
orders by the time-like Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution, 
while power terms of the form (rngZ/i^)" are negligibly small and can be safely neglected. In 
this massless- quark scheme oi zero-mass variable- flavour-number scheme (ZMVFNS), which 



is sometimes improperly referred to as NLL approximation j39|, the Q quark is treated as 
massless and appears as an active parton in the incoming hadron or photon, having a non- 
perturbative parton density function (PDF). The criterion ^ ^ mq is certainly satisfled for 
e"'"e~ annihilation on the Z-boson resonance, and for the photo-, lepto-, and hadropro duct ion 
of D and B hadrons with transverse momenta Pt ^ ttlq. Furthermore, the universality of the 
FFs is guaranteed by the factorization theorem , which entitles us to transfer information 
on how charm and bottom quarks hadronize to D and B hadrons, respectively, in a well- 
deflned quantitative way from e^e^ annihilation, where the measurements are usually most 
precise, to other kinds of experiments, such as photo-, lepto-, and hadroproduction. In 
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Refs. m 



12| |. the distributions in the scaled D- and S-hadron energy x = 2Ej^ measured 



at LEPl were fitted at LO and NLO in the ZMVFNS using, among others, the ansatz by 
Peterson et al. [15'] for the c ^ D and h ^ B FFs at the starting scale = 2mQ. In the 
D*^ {B+/B^) case, the e parameter was found to be Ec = 0.0851 and 0.116 [ej {sb = 0.0126 
and 0.0198 12]) at LO and NLO, respectively. We emphasize that the value of e carries 
no meaning by itself, but it depends on the underlying theory for the description of the 
fragmentation process, in particular, on the choice of the starting scale fiQ, on whether the 
analysis is performed in LO or NLO, and on how the final-state coUinear singularities are 
factorized in NLO. An alternative to the ZMVFNS with purely non-perturbative FFs is to 
decompose the FFs into a perturbative coirmonent, the so-called perturbative FFs (PFFs) 



decompose tne tts mto a perturbative coirmo 
[ifij ] , and a non-perturbative component [l^l • 



In the traditional massive-quark scheme or fixed- flavour-number scheme (FFNS), the Q 
quark is treated in the on-mass-shell renormalization scheme, as if it were a massive lepton 
in triplicate, and it only appears in the final state, but not as an active parton inside the 
incoming hadron or photon. There are no collinear singularities associated with the outgoing 
Q-quark lines that need to be subtracted and absorbed into FFs. This scheme breaks down 
for pt ^ niQ because of would-be collinear singularities of the form a^ln (p^^/j^g), which 
are not resummed. However, this scheme allows one to calculate a total cross section, which 
is infeasible in the ZMVFNS. Quantitative comparisons of the ZMVFNS and the FFNS for 
photoproduction in ep and 77 collisions may be found in Refs. Q and Q, respectively. 

A rigorous theoretical framework that retains the full finite-mg effects while preserv- 
ing the indispensible virtues of the factorization theorem, namely the universality and the 
DGLAP [13] scaling violations of the FFs, is provided by the general-mass variable-flavour- 
number scheme (GMVFNS) 2^. In a nutshell, this procedure consists in explicitly 
performing the rriq limit of the FFNS result, comparing the outcome, term by term, 
with the ZMVFNS result in the MS scheme, and subtracting the difference terms from 
the FFNS result. Owing to the factorization theorem |l^, the hard-scattering cross sec- 
tions thus obtained can then be convoluted with non-perturbative D- and i?-hadron FFs 
extracted from LEPl data using the pure MS scheme j6|,|l2|. This is consistent because the 
finite-mg terms omitted in Refs. ^, are relatively small, of order mg/m|. The impact of 
finite-mg terms on the proton PDFs was recently assessed by the CTEQ Collaboration {2^ . 
In this connection, we should also mention the so-called fixed-order next-to-leading-logarithm 
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(FONLL) scheme 



22j |. in which the ordinary result in the FFNS and a suitably subtracted 



22 



result in a ZMVFNS with PFFs are linearly combined using a certain weight function 
The conceptual merits of the GMVFNS and the FONLL are reviewed in Ref. j23| . 

The GMVFNS was recently implemented for direct 12411 a nd single-resolved j25| 77 colli- 
sions as well as for ep collisions with quasi-real photons [26[. In the case of 77 — D*^ + X 
at LEP2, the inclusion of finite-mc effects was found to reduce the cross section by approx- 



imately 20% (10%) at = 2m 



naively expected. From Refs. 



'3 mp) |24|. i.e., their magnitude is roughly m^/(p^)^, as 



25l |26| , we thus infer that finite-mc effects play a signifi- 
cant role only at rather small values, ~ 3 GeV, so that the ZMVFNS should yield a 
good approximation in the kinematic range > 3 GeV and > 6 GeV considered in a 
very recent ZEUS analysis 

The article is organized as follows. In Sec. II, we give a short description of the method. 
In Sec. Ill, we present a numerical analysis of -D*^ plus jet associated photoproduction, 



where we also compare to preliminary ZEUS data j9[. Section IV contains the conclusions. 



II. THEORETICAL FRAMEWORK 



Using the massless scheme, we can rely on the factorization theorem and write the pho- 
toproduction cross section for ep — > 7p D*^ + X as a convolution of the partonic cross 
section a with the PDFs of the incident particles and the FFs for an outgoing parton frag- 
menting into a D*^ meson. In this approach, the FFs are purely non-perturbative, universal, 
and subject to DGLAP [l^ evolution. We have 




X da'^-'^'^ixePe, XpPp, Pd/z, /i, M„ Mp, Mp), (1) 

where Mg and Mp are the initial-state factorization scales, Mp is the final-state factorization 
scale, and /x the renormalization scale. If a jet in addition to the D*^ meson is detected, a 
measurement function defining the jet has to be included in Eq. ([Q). The sum Ylijk I'uns 
over all partons, including quarks, gluons as well as photons, which can contribute if the 
energy of the subprocess is above their mass thresholds. Therefore, the charm quark also 
contributes as an incoming parton originating from the proton or the resolved photon. 
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At leading order, the subprocesses contributing to the partonic reaction ij kX can 
be divided into two categories, corresponding to a direct or a resolved photon in the initial 
state. The direct photon corresponds to i = 7 in Eq. (P), with F^/g approximated by the 
Weizsacker- Williams formula for the spectrum of the quasi-real photons, 

^ + a-y)\ ^ QL.a-y) 2(1 

y mly^ y \ ' 

where aem is Sommerfeld's fine-structure constant and Q^ax is the maximum virtuality of 
the photon. On the other hand, the resolved photon participates in the hard interaction 
via its quark and gluon content. In this case, Fi/e{xe, Me) is given by a convolution of the 
Weizsacker- Williams spectrum with the photon PDFs as 

Fi/e{xe,Me)= dydx^F^/e{y)Fi/^{x^,Me)S{yx^ - Xe). (3) 
Jo 

Typical examples of LO diagrams are depicted in Fig. ^ 




27r 



FIG. 1: Examples of LO diagrams involving (a) direct photons or (b) the quark or (c) gluon 
components of resolved photons. 

At NLO, an additional real or virtual parton can be radiated in the hard interaction. 
Representative examples of NLO diagrams are shown in Fig. |21 As the partons are massless, 
this leads to soft and collinear singularities. The soft singularities cancel between the real 
and the virtual corrections, while the collinear ones are absorbed into the PDFs or FFs. For 
example, the collinear singularities appearing at NLO if the incident photon splits into a 
collinear qq pair are absorbed into the PDF Fq/^{x^, M^) at the factorization scale M^. Thus 
the direct- and resolved-photon contributions separately exhibit strong dependences on Me, 
which cancel out only in their sum, up to terms which are formally beyond NLO. Therefore, 
it has to be stressed that, at NLO, only the sum of the two contributions carries a physical 
meaning. Technically, the infrared divergences have been isolated using a combination of 

5 



the phase-space-slicing [27'] and the subtraction [28[ methods. A more detailed description 
of the calculation can be found in Refs. Isoj] and will not be repeated here. The matrix 



elements 
EPHOX 



32 



lave been calculated in Refs. [31| and are implemented in the computer program 
which can serve to calculate the photoproduction of a large-p-r hadron or photon 
together with an optional jet. The program is constructed as a partonic event generator, so 
that the total cross section as well as differential distributions can be easily obtained. 



(a) 





(c) 



FIG. 2: Examples of NLO diagrams contributing to the real corrections in (a) direct and (b) 
resolved photoproduction, and (c) to the virtual corrections in direct photoproduction. 



III. NUMERICAL RESULTS 

We are now in a position to present our numerical analysis. We start by specifying our 
inputs and the kinematic situation. We work in the ZMVFNS with Uf = 5 massless quark 



flavors. For the PDFs of the proton, we take the MRST03 j33| set by Martin, Roberts, Stir- 
ling, and Thorne. For the PDFs of the photon, our default set is AFG04 13411 b v Aurenche, 
Fontannaz, and Guillet, which is an updated version of the original AFG parameteri- 
zation. In contrast to the AFG [3| set, the AFG04 jsJ] set also contains a bottom-quark 
PDF. The AFG04 Q PDFs are slightly higher at small values of x and lower at large 
values of x than the AFG [s^ PDFs, but the numerical difference is very small. In order to 
assess the potential of the ZEUS data 9] to constrain the photon PDFs, we also employ set 
GRV HO l36| by Gliick, Reya, and Vogt, which we transform from the DIS^ scheme to the 
MS scheme. For the D*^ FFs, we use the parameterization of Ref. p], where separate NLO 
fits to ALEPH and OPAL data are performed. As our default, we chose the set obtained 
from the fit to the OPAL data, as it has a lower value. However, the differences in 
the considered cross sections resulting from exchanging the two FF sets are negligible. For 



we use an exact solution of the two- loop renormalization group equation, where the 
asymptotic scale parameter A^^^ for rij = 5 is calculated from A*^^^ by requiring continuity of 
as at the bottom threshold /i = 2mfe, and A*^^^ is set to 278 MeV to be consistent with the 
MRST03 y proton PDFs. 



Our default scale choice is /i = mj' and M = 2mr, where mT = \/vnb\ + (pf?)^ is the 
transverse mass of the D*^ meson and we set vrtc = 1.5 GeV. As usual, we identify the 
three factorization scales Mg, Mp, and Mp in Eq. and denote their common value by 
M. In order to estimate the theoretical uncertainty in a conservative way, we vary fi and 
M independently about their default values, as /i = Cim^- and M = 2c2'mT with 1/2 < 
ci,C2 < 2. The maximum variation of the cross section thus obtained is shown as hatched 
green bands in Figs. OHSl For comparison, we also consider diagonal scale variations, where 
Ci = C2, in these figures. However, we caution the reader that the errors resulting from 
diagonal scale variations are too optimistic because of cancellations due to the fact that the 
dependences of the cross section on n and M act in opposite directions. 

We chose the kinematics in such a way that a direct comparison to preliminary ZEUS data 
on D*^ plus jet associated photoproduction is possible. These data have been produced 
with a proton energy of 920 GeV and an electron energy of 27.5 GeV in the laboratory frame, 
which corresponds to a ep center-of-mass (cm.) energy of y/s = 318 GeV. The maximum 
photon virtuality is Q^ax — 1 GeV^, and the photon-proton cm. energy W lies in the range 
130 GeV < W < 280 GeV. All rapidities rj refer to the laboratory frame, with the HERA 
convention that t he p roton is moving towards positive rapidity. The jets are defined using 
the fcy-algorithm 37[ with a jet radius of i? = 1. We employ the ZEUS convention 9] that 
an event is counted twice if two jets in addition to the D*^ meson are measured which both 
satisfy the cuts in transverse momentum and rapidity. However, we find that this case only 
occurs for about 0.3% of the generated events in the kinematic range considered. Unless 
stated otherwise, we adopt the acceptance cuts > 3 GeV, \ri^\ < 1.5, p:^ > 6 GeV, and 
—1.5 < rj^ < 2.4 from ZEUS 9]. We present histograms with the same binning as in Ref. jsj. 
As in Ref. P], we consider the sum of the cross sections for D*^ and D*^ mesons. 

We now present and discuss our figures. In Figs. OHil we confront the differential cross 

n 

section of ep D j + X in photoproduction as measured by ZEUS [9] with our NLO 
predictions. Specifically, Fig. El refers to da/dp-lp integrated over the full r]^ range (—1.5 < 
T]^ < 2.4) and Figs. IH Efa) and (b) to da/dr]^ integrated over the fuUp^ range (p^ > 6 GeV) 
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and the subintervals 6 < < 9 GeV and > 9 GeV, respectively. In addition, we 
present, in Figs. Eta)-(d), da/dplp integrated over the subintervals —1.5 < rj^ < —0.5, 
—0.5 < 7]^ < 0.5, 0.5 < 7]^ < 1.5, and 1.5 < r]^ < 2.4, for which experimental data are 
not yet available. In all cases, the D*^ meson is kinematically confined by the conditions 
p;^ > 3 GeV and \ri^\ < 1.5. The vertical bars on the ZEUS data |3] give the full errors, 
the purely statistical errors are indicated by the horizontal ticks on them. Our default 
predictions are represented by the solid red histograms and their errors due to independent 
scale variation by the hatched green bands. The errors due to diagonal scale variation 
are indicated by the dot-dashed and dotted pink histograms. The predictions evaluated for 
central scale choice with the GRV HO photon PDFs are given by the dashed blue histograms. 
We observe that our NLO predictions agree with the ZEUS data [9J within errors, except for 
two forward rj^ bins in the lower range, where the ZEUS data [9] slightly overshoot our 
NLO predictions. In the cases when the difference between the evaluations with the AFG04 
^1 and GRV HO photon PDFs is comparable to the experimental error, the GRV HO 
2S| set tends to yield a better description of the ZEUS data ^, except at large values of 
and rj^ . 

It is also instructive to look at the distribution da/dx]^^^, where the kinematic observable 



^obs is defined as 



X' 



p-f exp(— ?7 ) + p^ exp(— 77-'^ 



(4) 



obs — 2E'y 

As the contribution from the direct-photon subprocesses peaks at x^j^^ ~ 1, whereas resolved 
photons mainly contribute for x^j^^ < 1, a cut on x^j^g can serve to obtain samples enriched 
in direct- or resolved-photon processes. As already mentioned in Sec. II, the true direct- 
and resolved-photon contributions are related at NLO through factorization and, taken 
separately, exhibit strong Me dependences. Only their sum represents a physical observable 
that can be compared to experimental data. Notice that the individual parts can even be 
negative. In Fig. our central NLO prediction for the differential cross section da /dxl-^^ 
in the kinematic range p;^ > 3 GeV, \rj^\ < 1.5, p^ > 6 GeV, and —1.5 < rj^ < 2.4 (solid 
red histogram) is decomposed into its direct- (dashed blue histogram) and resolved-photon 
(dotted green histogram) components. We read off from Fig. [7| that, with our default scale 
choice, direct photoproduction dominates for x^-^^ > 0.7. 

In the case of single-hadron inclusive photoproduction at HERA, the direct- and resolved- 
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FIG. 3: The differential cross section da /dp!j, of D*^j+X in photoproduction for —1.5 < rf < 
2.4, > 3 GeV, and < 1.5 measured by ZEUS 9] is compared with om' NLO predictions (soUd 
red histogram) including conservative errors (hatched green bands). For comparison, the errors 
due to diagonal scale variation (dotted-dashed and dotted pink histograms) and the predictions 
evaluated for central scale choice with the GRV HO photon PDFs (dashed blue histogram) are also 
shown. 

photon contributions are known to be accumulated in the backward and forward directions, 
respectively .38]. It is interesting to find out if, in D*"^ plus jet associated photoproduc- 
tion, the rapidities t]^ and t]^ lend themselves as discriminators between direct and resolved 
photoproduction as well. To this end, we split our central NLO predictions for the differ- 
ential cross sections da/drj^ and da/dr]^ into their direct- and resolved-photon components 
and show the results in Figs. IHl^a) and (b) as the solid red and dashed blue histograms, 
respectively. The kinematic range considered for da/drj^ is p!j^ > 3 GeV, \ri^\ < 1.5, and 
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FIG. 4: Same as in Fig. |31 but for da/dij^ with > 6 GeV. 

plp> 6 GeV, while for da/drj^ it is > 3 GeV, > 6 GeV, and -1.5 < t]^ < 2.4. Notice 
that the superposition of the two histograms in Fig. |Hl^a) yields the solid histogram in Fig.^ 
We learn from Figs. Efa) and (b) that, in the kinematic range considered, the discriminat- 
ing power of the rapidity distribution with respect to direct and resolved photons, which is 
familiar from single-hadron inclusive photoproduction at HERA, carries over to da/drj^ , but 
not to da/drj^. 

We now explore the sensitivity of our NLO predictions of ep — > D*^j + X in photopro- 
duction to the gluon and charm-quark PDFs of the photon, so as to assess the potential of 
the HERA experiments to constrain these PDFs, which are presently less well known than 
those of the up, down, and strange quarks. In order to enhance this sensitivity, it is useful 
to suppress the direct-photon contribution. From the discussion of Figs. U\ and |Hl we know 
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FIG. 5: Same as in Fig. H but for (a) 6 < < 9 GeV and (b) > 9 GeV. 

that this can be achieved by focusing on low values values of x^^^^ and/or large values of 
rf ^ typically < 0.75 and rf < 0.5. Therefore, we reconsider in Fig. El the differential 
cross section da jdrf for > 3 GeV, < 1.5, and > 6 GeV (a) without and (b) with 
the acceptance cut x^y^^ < 0.75 (solid blue histograms) and turn off the gluon (dashed red 
histograms) and charm-quark (dotted green histograms) PDFs, one at a time. Notice that 
the solid histograms in Figs. El and IHl^a) are identical. We caution the reader that, strictly 
speaking, it is inconsistent to put to zero a photon PDF by hand because, in the determi- 
nation of the PDF set, it participated in the DGLAP |l^ evolution and thus influenced the 
other PDFs. Furthermore, its Me dependence is correlated with a similar Me dependence in 
the direct-photon contribution through the factorization procedure. Bearing these caveats 
in mind, it is nevertheless instructive to do so. We conclude from Fig. IHfa) and (b) that 
the contribution from the gluon inside the resolved photon is too small to be useful, while 
that from the charm quark is very significant. Quantitatively comparing Figs. EJ IH^a), and 
ini^a) after integration over i]^ , we find that the contribution due to the charm component in 
the photon makes up 92% of the resolved-photon contribution and 50% of the total cross 
section. After imposing the condition x]^^^ < 0.75, the latter fraction is increased to as much 
as 81%! In fact, the ZEUS data ^ in Fig. El overshoot the dotted green histogram in|Hfa) 
by several experimental standard deviations in the upper rj^ bins, which suggest that, in 
the framework of the ZMVFNS, the existence of intrinsic charm in the resolved photon is 
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FIG. 6: Same as in Fig. El but for (a) -1.5 < rf < -0.5, (b) -0.5 < rf < 0.5, (c) 0.5 < rf < 1.5, 
and (d) 1.5 < r]^ < 2.4 and without experimental data. 
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FIG. 7: Our central NLO prediction of the differential cross section da/dx'J^^^ of ep D*^j + X 
in photoproduction for pfj^ > 3 GeV, \r]^\ < 1.5, > 6 GeV, and —1.5 < r]^ < 2.4 (solid red 
histogram) and its direct- (dashed blue histogram) and resolved-photon (dotted green histogram) 
components. 

experimentally established. 

Another interesting observable is cos 6**, defined as 



As the angular dependence of subprocesses involving a gluon propagator in the t channel is 
approximately proportional to (1 — | cos^*|)~^, whereas it is proportional to (1 — | cos0*|)~^ in 
the case of a quark propagator, one can learn about the size of the contribution from diagrams 
of the type shown in Fig. ^b) by studying the differential cross section da /d cos 6*. A recent 
ZEUS analysis on dijet angular distributions in the photoproduction of open charm I7|] has 
shown that the measured cross section from resolved- enriched events, with x^^g < 0.75, 
exhibits a distinct asymmetry with a strong rise towards cos^^* = —1, i.e., the photon 
direction. This behaviour suggests that events with x^j^^ < 0.75 are dominantly produced by 
charm quarks coming from the photon side. On the other hand, the cos 6* distribution for 
direct- enriched events, with x]^^^ > 0.75, is almost symmetric, as expected for subprocesses 
like the one depicted in Fig. ^a). In order to substantiate these observations from the 
theoretical side, we now investigate the differential cross section da/d cos 9* of ep D*^i+X 




COS 9* = tanh 



(5) 
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FIG. 8: Direct- (solid red histograms) and resolved-photon (dashed blue histograms) components 
of our central NLO predictions of the differential cross sections (a) da/drj^ for > 3 GeV, 
\7]^\ < 1.5, and p^^ > 6 GeV and (b) da/dr]^ for > 3 GeV, p{, > 6 GeV, and -1.5 < ry-?' < 2.4 
of ep —>■ D*^j + X m. photoproduction. 

in photoproduction for > 3 GeV, \rj^\ < 1.5, > 6 GeV, and —1.5 < rj^ < 2.4 at NLO. 
In Fig. this cross section is decomposed in two ways: in the physically well-defined 
direct- (dotted blue histogram) and resolved-enriched (solid red histogram) contributions, 
with x^jjg > 0.75 and x^-^^ < 0.75, respectively; and in the mathematically defined direct- 
(dashed pink histogram) and resolved-photon (dot-dashed green histogram) contributions. 
From Fig. ^[ we observe that the direct- and resolved-enriched contributions, which can 
be measured experimentally, exhibit very similar cos9* dependences as their theoretical 
counterparts, which, taken separately, do not represent physical observables. In other words, 
the enriched contributions possess a rather high purity. Furthermore, we can confirm the 
findings of Ref. Q concerning the cos^* dependences of the direct- and resolved-enriched 
samples: the former is almost symmetric, whereas the latter is exhibits a steep rise towards 
the photon direction. This demonstrates that the bulk of the resolved-photon contribution 
is due to charm component, in accordance with our conclusions from Figs. Ufa) and (b). 
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FIG. 9: Our central NLO prediction of the differential cross section da /drj^ of ep D*^j + X 
in photoproduction for > 3 GeV, \r]^\ < 1.5, and pj, > 6 GeV (a) without and (b) with the 
acceptance cut x'^^^ < 0.75 (solid blue histograms) and the evaluations with Fg/^ = (dashed red 
histograms) or F^/^ = (dotted green histograms). 

IV. CONCLUSIONS 

Using results from Refs. we evaluated the cross section of D*^ plus jet asso- 

ciated photoproduction at HERA to NLO in the parton model of QCD implemented in 
the ZMVFNS with non-perturbative FFs extracted from LEPl lal, and studied various dis- 

Q 

tributions of it. As the most important result, we found that preliminary ZEUS data 
are nicely described by our theoretical predictions. This may be partly attributed to the 
fact that the acceptance cut p:^ > 6 GeV [9] excludes events at energy scales of order nic 
and below, which cannot be reliably described in this theoretical framework. In order to 
obtain a reliable prediction in the low-p^ range as well, the finite-m-c effects must be in- 
cluded, by working in the GMVFNS. In the case of D*"^ inclusive photoproduction, this 
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FIG. 10: Direct- (dashed pink histogram) and resolved-photon (dot-dashed green histogram) con- 
tributions as well as contributions with x^j^^ > 0.75 (dotted blue histogram) and x^j^g < 0.75 
(solid red histogram) to our central NLO prediction of the differential cross section da /d cos* of 
ep D*^j + X in photoproduction for > 3 GeV, \r]^\ < 1.5, pj, > 6 GeV, and —1.5 < r]^ < 2.4. 



was recently done for the direct-photon contribution in Ref. [26[. The good agreement with 

ri 

the ZEUS data [9] provides successful tests of the universality and the scalin g v iolations of 
the FFs, which are predicted by the factorization theorem and the DGLAP 13| evolution, 
respectively. 

Unfortunately, the variation of the NLO predictions due to scale changes is larger than the 
one stemming from using different contemporary NLO sets of photon PDFs, so that the latter 
cannot be further constrained by measurements of D*^ plus jet associated photoproduction 
at HERA. 

As for the relative importance of the various partons inside the resolved photon, charm 
was found to greatly dominate, while the gluon turned out to be practically irrelevant. In 
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particular, the dominance of the charm component in the resolved photon manifests itself in 
the characteristic cos 9* dependence of the cross section, a feature that was already exploited 
in experimental analyses j^l • The comparison of the ZEUS data [ol with our NLO predictions 
establishes the presence of intrinsic charm inside the resolved photon with overwhelming 
significance and thus supports the validity and usefulness of the ZMVFNS endowed with 
non-perturbative FFs in the kinematic regime considered. 



Acknowledgments 



We would like to thank E. Gallo, L. Gladilin, J.H. Loizides, and F. Sefkow for helpful 
discussions about Refs. and J.H. Loizides for making available to us the 

preliminary ZEUS data 9j in numerical form. B.A.K. thanks the Max-Planck-Institut fiir 
Physik for the hospitality extended to him during a visit when this work was finalized. This 
work was supported by the Bundesministerium fiir Bildung und Forschung through Grant 
No. 05 HT4GUA/4. 



[1] HI Collaboration, S. Aid et al, Nucl. Phys. B472, 32 (1996); HI Collaboration, C. Adloff et 

al, ibid. B545, 21 (1999). 
[2] ZEUS Collaboration, J. Breitweg et al, Phys. Lett. B 401, 192 (1997). 
[3] ZEUS Collaboration, J. Breitweg et al, Eur. Phys. J. C 6, 67 (1999). 

[4] S. Frixione, M.L. Mangano, P. Nason, and G. Ridolfi, Nucl. Phys. B412, 225 (1994); Phys. 

Lett. B 348, 633 (1995); S. Frixione, P. Nason, and G. Ridolfi, Nucl. Phys. B454, 3 (1995); 

S. Frixione and P. Nason, JHEP 0203, 053 (2002). 
[5] J. Binnewies, B.A. Kniehl, and G. Kramer, Z. Phys. C 76, 677 (1997); M. Cacciari and 

M. Greco, Phys. Rev. D 55, 7134 (1997); B.A. Kniehl, G. Kramer, and M. Spira, Z. Phys. C 

76, 689 (1997). 

[6] J. Binnewies, B.A. Kniehl, and G. Kramer, Phys. Rev. D 58, 014014 (1998). 

[7] ZEUS Collaboration, S. Chekanov et al, Phys. Lett. B 565, 87 (2003); L. Gladilin, in Proceed- 
ings of the XI International Workshop on Deep Inelastic Scattering (DIS 2003), St. Peters- 
burg, Russia, 2003, edited by Victor T. Kim (World Scientific, Singapore, 2004), Report No. 



17 



|hep- ex/0309043: S. Padhi, PhD Thesis, McGill University, Montreal, Canada, 2003, Report 

No. DESY-THESIS-2004-012. 
[8] ZEUS Cohaboration, T. Kohno for the cohaboration, in Proceedings of the XII International 

Workshop On Deep Inelastic Scattering (DIS 2004), Strbske Pleso, Slovakia, 2004 (to appear); 

T. Kohno, KEK-REPORT-2004-3, 2004. 
[9] ZEUS Collaboration, Abstracts 5-332 and 11-333, submitted to 32nd International Conference 

on High Energy Physics (ICIIEP'04), Beijing, China, 2004. 
[10] F. Sefkow, private communication; HI Collaboration, in preparation. 

[11] J. Binnewies, B.A. Kniehl, and G. Kramer, Z. Phys. C 65, 471 (1995); Phys. Rev. D 52, 4947 
(1995); B.A. Kniehl, G. Kramer, and B. Potter, Nucl. Phys. B582, 514 (2000); B597, 337 

(2001) ; S. Kretzer, Phys. Rev. D 62, 054001 (2000); L. Bourhis, M. Fontannaz, J.P. Guillet, 
and M. Werlen, Eur. Phys. J. C 19, 89 (2001). 

[12] J. Binnewies, B.A. Kniehl, and G. Kramer, Phys. Rev. D 58, 034016 (1998). 

[13] V.N. Gribov and L.N. Lipatov, Yad. Fiz. 15, 781 (1972) [Sov. J. Nucl. Phys. 15, 438 (1972)]; 

G. Altarelh and G. Parisi, Nucl. Phys. B126, 298 (1977); Yu.L. Dokshitser, Zh. Eksp. Teor. 

Fiz. 73, 1216 (1977) [Sov. Phys. JETP 46, 641 (1977)]. 
[14] J.C. Collins, Phys. Rev. D 58, 094002 (1998); private communication. 
[15] C. Peterson, D. Schlatter, I. Schmitt, and P.M. Zerwas, Phys. Rev. D 27, 105 (1983). 
[16] B. Mele and P. Nason, Phys. Lett. B 245, 635 (1990); Nucl. Phys. B361, 626 (1991); J.P. Ma, 

ibid. B506, 329 (1997); K. Melnikov and A. Mitov, Phys. Rev. D 70, 034027 (2004). 
[17] M. Cacciari, M. Greco, S. Rolh, and A. Tanzini, Phys. Rev. D 55, 2736 (1997). 
[18] B.A. Kniehl, M. Kramer, G. Kramer, and M. Spira, Phys. Lett. B 356, 539 (1995). 
[19] M. Cacciari, M. Greco, B.A. Kniehl, M. Kramer, G. Kramer, and M. Spira, Nucl. Phys. B466, 

173 (1996). 

[20] M.A.G. Aivazis, J.C. Collins, F.I. Olness, and W.-K. Tung, Phys. Rev. D 50, 3102 (1994); 

F.L Olness, R.J. Scalise, and W.-K. Tung, ibid. 59, 014506 (1999); A. Chuvakin, J. Smith, 

and W.L. van Neerven, ibid. 61, 096004 (2000). 
[21] S. Kretzer, H.L. Lai, F.I. Olness, and W.K. Tung, Phys. Rev. D 69, 114005 (2004). 
[22] M. Cacciari, M. Greco, and P. Nason, JHEP 9805, 007 (1998); M. Cacciari, S. Frixione, and 

P. Nason, ibid. 0103, 006 (2001); M. Cacciari and P. Nason, Phys. Rev. Lett. 89, 122003 

(2002) . 

18 



[23] B.A. Kniehl, in Proceedings of the l^th Topical Conference on Hadron Collider Physics: 
Hadron Collider Physics 2002, Karlsruhe, Germany, 2002, edited by M. Erdmann und Th. 
Miiller (Springer, Berlin, 2003), p. 161, Report No. hep-ph/0211008 ; B.A. Kniehl and F. Se- 
fkow, in Proceedings of the XI International Workshop on Deep Inelastic Scattering (DIS 

2003) , St. Petersburg, Russia, 2003, edited by Victor T. Kim (World Scientific, Singapore, 

2004) , Report No. |hep-ph/03 12054] 

[24] G. Kramer and H. Spiesberger, Eur. Phys. J. C 22, 289 (2001). 
[25] G. Kramer and H. Spiesberger, Eur. Phys. J. C 28, 495 (2003). 



[26] G. Kramer and H. Spiesberger, Report No. MZ-TH/03-18 and |hep-ph/03 1 1062| (unpubhshed) . 

[27] W.T. Giele and E.W.N. Glover, Phys. Rev. D 46, 1980 (1992); W.T. Giele, E.W.N. Glover, 
and D.A. Kosower, Nucl. Phys. B403, 633 (1993). 

[28] S. Frixione, Z. Kunszt, and A. Signer, Nucl. Phys. B467, 399 (1996); S. Catani and M.H. Sey- 
mour, ibid. B485, 291 (1997); B510, 503(E) (1997). 

[29] P. Chiappetta, R. Fergani, and J.P. GuiUet, Z. Phys. C 69, 443 (1996); M. Fontannaz, 
J. P. Guillet, and G. Heinrich, Eur. Phys. J. C 21, 303 (2001); 22, 303 (2001); M. Fontannaz 
and G. Heinrich, ibid. 34, 191 (2004). 

[30] M. Fontannaz, J.P. Guillet, and G. Heinrich, Eur. Phys. J. C 26, 209 (2002). 

[31] P. Aurenche, A. Douiri, R. Baier, M. Fontannaz, and D. Schiff, Z. Phys. C 24, 309 (1984); 
R.K. Ellis and J.C. Sexton, Nucl. Phys. B269, 445 (1986); P. Aurenche, R. Baier, A. Douiri, 
M. Fontannaz, and D. SchifF, ibid. B286, 553 (1987). 

[32] M. Fontannaz, J.P. Guillet, and G. Heinrich, the computer program can be downloaded from 
the URL |http : //wwwlapp . iii2p3 . f r/lap th/PHDX JAMILY/ main . ht ml see also Ref. j3o|. 

[33] A.D. Martin, R.G. Roberts, W.J. Stirling, and R.S. Thorne, in Proceedings of the XI Interna- 
tional Workshop on Deep Inelastic Scattering (DIS 2003), St. Petersburg, Russia, 2003, edited 
by Victor T. Kim (World Scientific, Singapore, 2004), Report No. IPPP/03/43, DCPT/03/86, 



Cavendish-HEP-2003/13, and |hep-ph/03 07262 , 
[34] P. Aurenche, J.P. Guillet, and M. Fontannaz, private communication. 
[35] P. Aurenche, J.P. Guillet, and M. Fontannaz, Z. Phys. C 64, 621 (1994). 
[36] M. Gliick, E. Reya, and A. Vogt, Phys. Rev. D 45, 3986 (1992); 46, 1973 (1992). 
[37] S. Catani, Yu.L. Dokshitzer, M.H. Seymour, and B.R. Webber, Nucl. Phys. B406, 187 (1993); 

S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993). 



19 



[38] F.M. Borzumati, B.A. Kniehl, and G. Kramer, Z. Phys. C 59, 341 (1993); B.A. Kniehl and 

G. Kramer, ibid. 62, 53 (1994). 
[39] The non-logarithmic corrections of relative order ag are fully included, except for terms that 

are suppressed by powers of frng///^^ . 



20 



